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Osl ■ Abstract 

S ■ A notion of stochastic deformation is introduced and the corresponding algebraic deformation procedure 

is developed. This procedure is analogous to the deformation of an algebra of observables like deformation 
^^ ' quantization, but for an imaginary deformation parameter (the Planck constant). This method is demon- 

strated on diverse relativistic and nonrelativistic models with finite and infinite degrees of freedom. It is 

r^ ■ shown that under stochastic deformation the model of a nonrelativistic particle interacting with the elec- 

o ■ 

q_) ' tromagnetic field on a curved background passes into the stochastic model described by the Fokker-Planck 

H ! equation with the diffusion tensor being the inverse metric tensor. The first stochastic correction to the 

j^ , Newton equations for this system is found. The Klein-Kramers equation is also derived as the stochastic de- 

t/^ . formation of a certain classical model. Relativistic generalizations of the Fokker-Planck and Klein-Kramers 

C^ ' equations are obtained by applying the procedure of stochastic deformation to appropriate relativistic classi- 
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cal models. The analog of the Fokker-Planck equation associated with the stochastic Lorentz-Dirac equation 
is derived too. The stochastic deformation of the models of a free scalar field and an electromagnetic field 



O i is investigated. It turns out that in the latter case the obtained stochastic model describes a fluctuating 



electromagnetic field in a transparent medium. 
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I. INTRODUCTION 

Stochastic equations are usually considered as an effective tool that is able to take into account 
infinitely many factors acting on an open system. Generally, they are obtained by adding to the 
classical equations of motion of the system a random force (noise) with some probability distribution 
law. In spite of this common view we shall develop here a slightly different technique inspired by 
an algebraic approach to quantum mechanics. As is known quantum fluctuations or quantum 
mechanics itself can be regarded as resulting from a deformation of underlying classical mechanics 
or, more precisely, from an algebraic deformation of the Poisson structure on the phase space of the 
classical system. This approach, known as deformation quantization, was initiated in the seminal 
papers |Q]. Since the quantum fluctuations arise from the deformation of the Poisson structure, 
it is reasonable to pose the question whether the stochastic fluctuations could be reproduced in 
such a way. In this paper we give an affirmative answer on this question and sustain it by diverse 
examples. 

A similarity between equations of quantum and stochastic mechanics was perceived by many 
authors starting from Schrodinger himself Q. Much afterwards this similarity was brought in an 
almost perfect accordance by Zambrini in his works [^ on Euclidean quantum mechanics^. However 
the key point of these works as well as the works of other authors [§, |5|, y, 0, ^ was an analytic 
continuation to imaginary time and not the deformation of the Poisson structure. Further we shall 
see that in the relativistic case this slight disagreement is not a matter of interpretation but gives 
inequivalent results. 

The principal point of this paper is that we can obtain stochastic mechanics deforming the 
associative and commutative product of the algebra of classical observables (the smooth functions 
over the phase space) along the same lines as in the deformation quantization procedure. However 
contrary to an ordinary quantum mechanics with the real deformation parameter (the Planck 
constant) we deform the product by an imaginary deformation parameter. By analogy we shall 
refer to this procedure as stochastic deformation^. The deformation parameter, which we denote by 
h, characterizes the "openness" of the system and a variance of random forces is proportional to it. 
The mention should be made that the Hamiltonian generating the dynamics of stochastic system 



^ It is noteworthy to say that unlike Zambrini's works we do not construct the resulting stochastic mechanics as a 
solution of the so-called Schrodinger gedanken experiment but expound it in a rather traditional way. 
We do not use the more familiar term "stochastic quantization" by two reasons. Firstly, this term refers to the 
fixed notion of theoretical physics which is not concerned with what we are considering in this paper. Secondly, 
there is nothing quantal in the theory developed here. That is to say such a deformation does not result in a 
quantization of observables. 



may take forms unexpected at first glance, because the classical limit ^ ^ of stochastic mechanics 
is not an ordinary classical mechanics. In this limit we obtain some classical mechanics where the 
momenta should be regarded as systematic forces (or proportional to them) acting on the system, 
for this limit is equivalent to a treatment of the system in the equilibrium state when stochastic 
fluctuations vanish and the effect of systematic forces is balanced by a dissipation. The expansion 
of the pq-syrahol of the Hamiltonian in momenta is closely related to the Kramers-Moyal expansion 
and the coefficients of this expansion are connected with the cumulants of the probability density 
function of a noise. In the simplest case of a linear symplectic space the quadratic in momenta 
Hamiltonian corresponds to a Gaussian noise. 

To avoid misunderstanding we note in advance that the result of stochastic deformation of a 
classical model, that we call stochastic mechanics, is not Nelson's stochastic mechanics Q developed 
to give a stochastic interpretation of quantum mechanics. As a practical matter, our approach can 
be regarded as a gauged version of the operator approach to stochastic mechanics expounded in |^, 
g, 0, P, IC, ^. Such a formulation offers the advantage that one can apply the developed methods 



of quantum mechanics and field theory to stochastic mechanics almost without any changing. In 
this context it is sufficient to mention, for instance, the models with gauge symmetries. The 



methods of quantization of systems with constraints (see for the introduction |12, 13, |J]) adapted 
to stochastic deformation give the technique of generation of stochastic models respecting not 
merely the global symmetries of the initial classical system but the gauge symmetries as well. In 
other words, these methods allow us to introduce the noise to physical degrees of freedom that are 
not expressed in an explicit form. We shall consider several models with gauge symmetries in this 
article. 

The paper is organized as follows. We start our investigation with a formulation of general rules 
of stochastic deformation by the example of a linear symplectic space (Section ||) . In developing this 
scheme we follow the central proposal of deformation quantization and try to devise the procedure 
in a more algebraic way. We do not deal with the linear space of states and its dual but with the 
algebra of operators on them which are recognized as stochastic observables. A state of a stochastic 
system is specified by an analog of the density operator. For pure states it is useful to realize the 
algebra of observables in some linear space, at that we construct this space by the action of creation 
operators on the vacuum state. This standpoint allows us to calculate averages and make some 
proofs by means of the basic relation of the Heisenberg-Weyl algebra only. Generally speaking, the 
realization of some operators arising on intermediate steps of calculations can lead to divergent 
integrals or series that cancel each other in the final result. Then the algebraic approach can be 



considered as some regularization or prescription to handle these singular integrals. Besides, on 
this algebraic way we can establish stochastic mechanics for systems with nonlinear phase space 
such as symplectic or even non-regular Poisson manifolds. Explicit constructions of the deformed 



product (star-product) of the algebra of smooth functions on such manifolds are given in |15, 16|. 
In conclusion of Section || we relate the mechanics obtained from stochastic deformation and 
ordinarily formulated stochastic mechanics using a path-integral representation of the transition 
probability and the Langevin equations associated with it. There we present the formal relation 
only. The existence problems are left beyond the scope of the paper. The very method of stochastic 
deformation implies that the partial differential equations of the Fokker-Planck type are primary 
for this framework, while the stochastic equations of the Langevin type are only used to give 
a more lucid physical interpretation to the obtained equations. At the same time the form of 
the Hamiltonian generating an evolution is related to the probability density function of the noise. 
Thereby the information about the probability distribution law for the noise enters to the stochastic 
deformation procedure. 



In Section III we consider several significant models revealing the key features of the developed 



formalism. All of them are divided into three categories. The first class exposed in Section HI A 
is constituted by the models of a nonrelativistic particle under the influence of a random force. It 
is a classical subject of stochastic mechanics and we consider these models to include them into a 
general scheme and gain an experience required in understanding of subsequent sections. Namely, in 
this section we study the stochastic deformation of the models of a nonrelativistic particle coupled 
to the electromagnetic and gravitational fields as well as the stochastic deformation of the model 
leading to the Klein-Kramers equation. Under stochastic deformation the first two models result 
in the stochastic systems described by the Fokker-Planck equations with trivial and nontrivial 
diffusion tensors respectively, the diffusion tensor being the inverse metric [o]. 



In Section IIIB we extend our analysis to relativistic models with constraints. We investigate 
three stochastic models of a relativistic particle affected by a random force. The first model is 
the stochastically deformed model of a relativistic particle interacting with the electromagnetic 
field. We show that such a model is described by a relativistic equation which generalizes the 
Fokker-Planck equation in the same sense as the Klein-Gordon equation generalizes the Schrodinger 
equation. For a constant systematic force acting on the particle the obtained equation is in one- 
to-one correspondence with the relativistic diffusion equation ||l^. Making use of the descent 
method [l^, |2^ we derive a path-integral representation of the transition probability, whence it 
turns out that such a relativistic Fokker-Planck equation corresponds to a relativistic particle 



under the influence of a non-Gaussian noise. The second model examined in this section is a 
relativistic generahzation of the model leading under stochastic deformation to the Klein-Kramers 
equation. The respective relativistic equation proves to describe a relativistic diffusion studied 



in [21, ^, ^, 24, 25, 26]. Here we also obtain a path-integral representation of the transition 
probability both in the gauge of a laboratory time and in the proper time gauge. The third 
model of this section is related to stochastic mechanics of a relativistic charged particle with 
the radiation reaction taken into account, i.e., we derive in this section a relativistic analog of 
the Fokker-Planck equation associated with the stochastic Lorentz-Dirac equation. It looks very 
likely that the procedure developed in these sections is straightforwardly generalized to stochastic 
reparameterization invariant relativistic ordinary differential equations of an arbitrary order. 



Section IIIC is devoted to stochastic deformation of relativistic field theories. We only touch the 
problem and consider free models of the scalar and electromagnetic fields. Under the assumptions 
of causality and relativistic invariance we derive the propagators for these models. In the case of 
the electromagnetic field the obtained propagator coincides with the well-known correlator of the 



electromagnetic fields in a transparent medium |27|. Besides the last model illustrates how the 



BRST-quantization technique can be applied to stochastic mechanics. 

Thus we see that the stochastic deformation method covers a considerable part of stochastic 
physics and apparently any other stochastic model can be formulated in terms of this unifying 



algebraic approach. The first model in Section IIIB shows that the non-Markovian processes can 



be obtained by stochastic deformation as well. The non-Markovian processes have to be described 
by the models with constraints. Inasmuch as we regard in this paper many topics of theoretical 
physics the reference list is, of course, incomplete. I tried to make references to the basic works 
known to me and to the works that might be useful to the reader in understanding and further 
developing of the scheme evolved here. 

We use the following notation and conventions. The system of units is chosen so that the 
velocity of light c = 1 . Greek letters denote space-time indices and Latin indices indicate the spatial 
components of tensors. Sometimes we shall use boldface characters to denote the spatial part of 
coordinates. Einstein's summation rule is assumed unless otherwise stated. Usually overdots will 
denote a differentiation with respect to time. In the sections regarding nonrelativistic models d is 
a dimension of the configuration space with the coordinates x and indices are risen and lowered by 
the metric tensor ry^i^ = diag{—l, 1, . . . , 1), while for the relativistic models d is a dimension of the 
space-time, x is the space-time coordinates and the metric tensor rj^i, reads as diag{l, —1, . . . , —1). 



II. THE RULES OF STOCHASTIC DEFORMATION 

In this section we formulate the rules of stochastic deformation and establish the relation be- 
tween a stochastically deformed model and a certain stochastic mechanics'^. 

Consider a classical system on the linear symplectic manifold M with canonical coordinates x* 
and pj : 

{x\pj] = 5i, i,j=T^, (1) 

where d is a dimension of the configuration space and curly brackets denote the Poisson brackets. 
The algebra of classical observables is the real commutative associative algebra of smooth functions 
on M. The evolution is generated by the observable H{t,q,p) that is the Hamilton function. 

Then we deform the algebra of classical observables in the manner of deformation quantization 
but with an imaginary deformation parameter (the Planck constant) such that 

[x\pj] = hsi, h>o. (2) 

Hereinafter we denote elements of the deformed algebra by hats and imply the Weyl-Moyal star- 
product^ 

/^ = E ;^ 2 ''"''" • • • ^""'"^'^--n/(^)56,....„5(^), (3) 

n=0 ■ ^ ^ 



where z = {x,p), an,bn = I, 2d, the functions /(z) and g{z) are the Weyl symbols of the corre- 
sponding elements of the deformed algebra, u> is the inverse to the symplectic 2-form cuab- Recall 
that the Weyl-ordered operator for a monomial in momenta symbol looks like 

1 " 

V'^-'-{x)pi,...pi„ ^ —Y,CnPil■■■P^,V"'■■""{x)p^,+,■..p^r., (4) 

fe=0 

where C^ are the binomial coefficients. The physical meaning of the constant h in formula (|2|) will 
be elucidated below. Roughly, h characterizes a variance of stochastic forces acting on a classical 
system. 

Let us given a linear functional Sp on the deformed algebra mapping to real numbers and 
vanishing on commutators, viz. 

Sp(/5) = Sp(5/), V/,5, (5) 



The general scheme evolved in this paper was sketched in the preprint [ESf. 
* For the introduction to different orderings, symbols and star-products see, e.g., 



which we shall call the trace. An explicit formula for the trace of the element / has the form 






Sp/= / TTTZ^fi^'ip) ^ Mf9)= / j:-T;rjfix,ip)g{x,ip). (6) 



Then we define a complete set of elements {p\} of the deformed algebra by the properties 

A 

2. Sppx = l, ^^^ 

3. pxpy = Sxx'px- 

Consider a class of the complete sets related to each other by similarity transformations, i.e., 
two complete sets {px} and {o"^} are in the same class iff there exists an invertible element U in 
the algebra spanned on the generators x* and pj such that 

^/(A) = U-'pxU, (8) 

where / is a bijection and U does not depend on A. We choose the class which contains the 
complete set {px} corresponding to the x-representation 



xdf-i i\ I d p 

Px = 6 {x -x)= I 7:^3-^ exp 



^Pi{x' - x') 



(9) 



{27rh) 

This class also contains the complete sets associated with any Lagrangian section of the symplectic 
space M obtained from the coordinate Lagrangian section by a linear symplectic transformation. 
The element U realizing a similarity transformation is a solution of the equation 

hU{t) = ^z'^iOabAlz'^Uit), U{0) = i, (10) 

for appropriate t, where the matrix A^ belongs to the Lie algebra of the symplectic group. 

We say that the element T of the deformed algebra is a stochastic observable iff there exists a 
complete set {tx} from the chosen class such that 



f = Y,nx)ix, (11) 



A 

where r(A) is a certain real-valued function. In other words the stochastic observables should be 
diagonalizable in the chosen class. 

The state of the stochastic system is characterized by the observable p with a unit trace: 

Spp = l. (12) 



The pure state is specified by an additional idempotency requirement 

p2 = p. (13) 

The average of an observable T over a state p is defined as 

(f) := Sp(/5f ). (14) 

In particular, the average over the state /5 of a certain pure state p\ gives the probability to find a 
stochastic system in this pure state. 

The dynamics of a stochastic system in the state p are generated by the Weyl-ordered operator 
H corresponding to the Hamilton function H{t,x,p) and described by the von Neumann equation 

h/)=[H,p], (15) 

whence the "Euclidean" Heisenberg equation for averages follows 

h^^{f) = {hdtf+[f,H]). (16) 



The evolution defined in this way maps observables into observables. Besides it follows from (15) 
that the evolution of a probability density function to find a system in the pure state px obeys the 
equation 

hj^{px)p = Sp{px[H,p]), (17) 

which is nothing but the Fokker-Planck equation. The conservation of the total probability is a 
consequence of the trace property. 

Now we are in position to touch the problem of an arbitrariness in a definition of the star- 
product related to the ordering prescription. This arbitrariness at least for a linear symplectic 
space does not affect averages and their evolution, but as in quantum mechanics when we convert 
one star-product to another some observables becomes nondiagonalizable in a class of complete 
sets corresponding to the new star-product. Anyway different orderings for the same physical 
observable or the Hamiltonian just result in different definitions of the quantities measured in 
stochastic mechanics (correlators) in terms of the coefficients of an expansion of the observable and 
the Hamiltonian in momenta. 

Since the phase space of the given classical system is a linear symplectic space it is useful to 
realize the deformed algebra by operators acting in some linear subspace V of the linear space of 



smooth functions on the configuration space. We introduce the Dirac notations for elements from 
the hnear space V and its duah 

IV) e V, {ifl G V*, (18) 

where V* is the hnear space of hnear functionals on V. Operators acting in the space V are 
naturahy translated to operators acting in the dual space V*. The eigenvectors from the space V 
and V* of the operator T corresponding to the eigenvalue t we denote by 

IT = t), or It), and (T = t\, or {t\, 

(19) 

T\t)=t\t), {t\T={t\t. 

For the given left and right vectors (vacua) \x = 0) and {x = 0| we construct the eigenvectors 
of the operator x 

\x = a) = eii'''P'\x = 0), {x = a\ = {x = 0\e-T.''*P\ (20) 

Their inner product possesses the properties 

f{x){x\x') = f{x'){x\x'), yfeV, (x + A|x' + A) = (x|x'), VA, (21) 

and, consequently, is equal up to a constant factor to the (5-function. On setting this constant to 
unity the element of the complete set @ can be written as 

Px = \x){x\. (22) 

In the same manner we can prove that an element of a complete set associated with any other 
Lagrangian section has the from (p2|) with appropriate eigenvectors. The linear space V is spanned 
on the vectors obtained from \x) by an action of the various operators U realizing the similarity 
transformation in the chosen class of complete sets. 
Thus for the pure state 

p=\ij){^\, (^|V) = 1, (23) 

the probability to find a system in the point x at the time t takes the form 

p{t,x) = {x\^){^\x). (24) 



Now it is easy to obtain a formal solution to the Fokker-Planck equation (^). Its fundamental 
solution or, in physical terms, the transition probability can be represented in two ways^ 

Git', x'; t, x) := * ' V ' = '' \ ' '' \ (25) 

where Ut' ^t is the evolution operator obeying the equations 

ndt,Ut',t = H{t')Ut',u Ut,t = i. (26) 



A convolution of ( p5[ ) with (|2^) gives the probability density function at the moment t' . To provide 
linearity of equations ( [T7| ) with respect to p{t, x) we have to claim that one of the kernels ( p5[ ) is 
independent of p(i,x). Therefore either {ip\x) or {x\ip) is independent of p{t,x). 
Consequently, the density operator (|2^) in the x-representation looks like 



{x\i;){^\y) = ei[^(*'^)-^(*'-)]p(t,x), (xmiy^ly) = ei[^(*'^)-^(*'-)lp(t,y), (27) 

for the first and second cases respectively. Here we introduce S{t, x) which is the analog of a 
quantum mechanical phase ||3|, ^]. This function can have discontinuities or even be complex but its 
gradient contributing to the observable averages should be real and have removable discontinuities 
only. The first representation in ( p5| ) for the fundamental solution is a forward evolution operator 
twisted by e^^ '*'^^ and it takes the form of a forward transition probability, while the second 
representation of the fundamental solution has the form of a backward transition probability for 
the Fokker-Planck equation. 

These transition probabilities depend on the phase having different definitions for the first and 
second cases (p7|). For the given initial probability density function one can choose the initial 
phase in such a way that the forward transition probability will give the same probabilities as the 
process generated by the backward transition probability with the same initial probability density 
function. The phase is a new "degree of freedom" as against to classical mechanics and the above 
mentioned flexibility of the formalism is only concerned with its interpretation in terms of stochastic 
mechanical quantities. For the forward transition probability with the definition of phase given 
by the first formula in (|27| ) this interpretation takes a simplest form and will be discussed below. 
Henceforth we shall refer to the forward transition probability as the transition probability. 



^ Other representations are permissible if one factorizes the 5- function, i.e., 

S{x~y) = if{x,y)ip{x,y), 

where ip{x,y) and tp(x,y) are not proportional to the (5-function. This is possible, for example, in the multidimen- 
sional case or by the use of exotic functions such as the fractional powers of the 5-function (see, e.g., ISOl). In both 
cases one should take care about linearity of Eq. ( |l7[ ) with respect to p{t,x). 

10 



A. Stochastic deformation and Langevin equation 

By standard means (see, e.g., [p9| , [84 ) we can construct a path- integral representation of the 
transition probabihty ( p5[ ) . Notice in advance that we shall not consider here the existence problems 
of fundamental solutions and their well-defined path-integral representations. The interested reader 
can consult the references above. 

Using the relation between the kernel of the Weyl-ordered operator T in the coordinate repre- 
sentation and its symbol 

we formally have 

{^{t + dt)\x'){x'\Ut+dt,tW 



{^{t)\x) 
(x'lexpi— H{t,x,p + VS{t,x))+dtS{t,x) \ \x) = 

^'P(*) exp |-^ [pi{t)±\t) + i{Hw {t,x{t),w{t)) + dtS{t,3:{t)))] dt }, (29) 



(271^)^^ ^1 h 

where Hwit,x,p) is a Weyl-symbol of the Hamiltonian H in which the momenta operators pj are 
replaced by pj + djS and 

x{t) = x, x{t + dt) = x', x{t) := {x(t + dt)-x{t))/dt, x{t) := {x{t + dt) + x{t))/2. 



As long as the transition probability ( |25D possesses the defining property of the Markov process we 
can cut the time interval \t, t'] into pieces for which formula (^) makes sense and then integrate 
over intermediate positions. As a result the fundamental solution to the Fokker-Planck equation 
( [ITD takes the form 

d'^pir) 



G{t',x';t,x)=J n dMr) fl 



(27rh¥ 

Te{t,t') rG[t,t') ^ ^ 

expl-- J dT[pi{T)±'{T)+i{Hw{T,x{T),ip{T))+drS{T,x{T)))]}. (30) 

Notice that the main contribution to the transition probability (pH) is made by paths approximating 
a classical trajectory. This becomes manifest if one makes the change of variables pj — > pj — djS 
and neglects stochastic corrections. 

On integrating (pO) over momenta we arrive at the functional integral with the expression in 
the exponent taking the form of the Stratonovich-type stochastic integral (see, e.g. [^). It is a 

11 



consequence of the use of Weyl-symbol for the evolution operator. If we had chosen the pg-symbol, 
we would have obtained the Ito-type stochastic integral. Namely, in that case the transition 



UX 



probability (30) becomes 

Git',.';t,x)= f n dMr) n S£^ 

rg(t,t') re[t,i') ^ ' 

f ^ 't 

expl-- J dT[pi{T)±\T)+i{H{T,x{T),ip{T))+drS{T,x{T)))]\, (31) 

where H(t, x,p) is a pg-symbol of the Hamiltonian operator with the momenta pj + djS. Certainly, 
the transition probabilities ( pOD and (31) are equal to each other as they are different symbols of 
the same evolution operator. 

Let us now provide an interpretation of the above stochastic mechanics in terms of the Langevin 
equation 

x\t) - v\t) = 0. (32) 

To this end we insert the 5-function with the LHS of (p^) in its argument into the path-integral 



(31) and integrate it over i^*(t). Then 

F{v,T,x) := JJ^e-^'^^P^'^'exp!^^ [F(r,x,ip) + a,5(r,x)] | , (33) 

can be interpreted as a probability density function of the noise i^*(t) at the time slice r, the 
Langevin equation ( |32D being understood in the Ito sense. The function 

$(A, r, x) := / d'^ue'^^'''F{u, r, x) = {dr)-'^ ^^P | ^ ^ (^' ^' '^7") + ^^^^'^' ^^ I ' ^^^^ 

is known as the characteristic function and the coefficients of the Taylor series in iX of the expression 
in the exponent are called the cumulants. For example, the second cumulant is proportional to the 
inverse metric (inverse mass matrix) in the Hamiltonian H{t,x,p) multiplied by the deformation 
parameter h. The second cumulant is equal to the mean-squared deviation that is why we claimed 
that h characterizes the variance of stochastic forces acting on the system. 

Now it is easy to see that the Hamiltonian at most quadratic in momenta leads to the Langevin 
equation with a noise whose the probability density function is a product of delta and Gaus- 
sian functions. The Hamiltonians depending on higher powers of momenta correspond to a non- 
Gaussian noise. If the Hamiltonian is analytic in momenta then in the classical limit, fi — > 0, only 
the first cumulant survives and the transition probability ( |3ll ) becomes Liouvillian. The zeroth 

12 



cumulant is a mere normalization factor and, as follows from the normalization condition on (31), 
is equal to —din dr. 

The above interpretation in terms of the Langevin equation works well when the probability 
density ( |33|) of noise is a positive and normalizable function. If this function takes negative values 
then the transition probability (^) also possesses negative values for some values of its arguments. 
However it does not automatically imply that the considered stochastic mechanics is unphysical. 
For example, this fact may signify that the (5-localized probability density functions are not well 
defined for such a system from the physical point of view. At the same time there may exist a 
class of probability density functions for which a convolution with the transition probability gives 
reasonable results. We shall encounter with the problem of negative probabilities in considering a 



relativistic generalization of the Fokker-Planck equation in Section III B l|. 



III. EXAMPLES 

A. Nonrelativistic particle 

In this subsection we consider the stochastic deformation of three nonrelativistic models: a 
particle interacting with the electromagnetic field, the same model on a curved background and 
the model leading to the Klein-Kramers equation. For these models we obtain a path-integral 
representations of the transition probabilities and the associated Langevin equations. Some simple 
applications of the developed scheme are also demonstrated. 

1. Fokker-Planck equation 

According to general rules expounded in the previous section the Weyl-ordered Hamiltonian for 
a nonrelativistic particle looks like 

H=^ll^^ + A''. (35) 

It is convenient to realize the operators x* and pj in the linear space V of functions on the config- 
uration space in the following way 

x' = x\ pi = -hdi. (36) 



13 



Then for the pure state of the form (p3| ) the von Neumann equation (p^ ) reduces to the two 
"EucUdean" Schrodinger equations 



hdtilj{t,x) 



jPi - Ajf 
2m 



An 



i){t,x), hdtip{t,x) 



{Pi + Aj 
2m 



An 



ipit,x), (37) 



provided that the standard inner product is understood. For a correct stochastic interpretation 
the functions ip{t,x) and ip{t,x) have to be both positive (or both negative). The fields A^ are the 
gauge fields, which we shall call the electromagnetic fields. Their physical meaning is obvious from 
the general considerations of the previous section. Namely, introducing the phase, S := h\inp, we 
see from (b3) that the quantity 



diS - Ai 



(38) 



is proportional to the first cumulant of the probability density function of noise ( [3^ ) whereas the 
first cumulant is equal to the expectation value of a random variable. We dwell on the interpretation 



of A^ a bit later. 



The Schrodinger equations (37) are invariant under the following gauge transformations 
V'(t,x)^V(i,a:)e-«(*'^), (/j(i,x)^(^(t,x)e«(*'^), A^{t,x) ^ A^{t,x) + df,C{t,x). (39) 

In particular, these transformations do not change the probability density function. The conserved 



4-current corresponding to the gauge transformations (39) is given by 



j^^ = ( (/.V, ^ W - A')^ - ^l^if + A')ip\ ] . 



(40) 



The system of equations ( |37D is Lagrangian with the Hamiltonian action of the form 



SH[^,ij] = j dt{if\hdt - H\ 



(41) 



that is the fields ^{t,x) and ip{t,x) are canonically conjugate with respect to the Poisson bracket. 



On substituting the definition of phase (27), the system of evolutionary equations ( |37|) can be 
rewritten in the equivalent form 



dtp = -d' 



h d,s - A, ■ 

-7^dip+ p 

2m. m 



dtS -Ao + 



^^'^ ^'^' ^d\d.S-A,). (42) 



2m 



2m, 



The first equation is Fokker-Planck equation, which is the divergence of the 4-current (|40|), while 
the second equation can be referred to as the quantum Hamilton- Jacobi equation [pi] , |3^ . Now it 
is evident that if one neglects stochastic corrections then the initially (5-shaped probability density 
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function p{t, x) keeps its own form and propagates as a classical charged particle in the electro- 
magnetic field^ with particle's momentum diS{t,x) — Ai{t,x). 

Taking into account the stochastic corrections we see that the more the probability density 
function is localized, the higher is the probability flow resisting localization. It can be perceived, 
for example, from the stochastic analog of the quantum mechanical uncertainty relation 

{{xr)Msf) > Y' (43) 

where summation is not understood, p^^ := —hd^lnp^'"^ is the osmotic momentum. This relation 
is easily deduced from the inequality 

/ d'^x U^x' - na,)p^/^j ^ > 0, VC G M, (44) 

under the assumption that p{x) tends to zero at spatial infinity faster than x^^. 

To find the first stochastic correction to the classical equations of motion we use the Heisenberg 



equations ([iq) for averages of the operators x* and pj 



m^ix,) = {p,-Ai) = {d,S~A,), ^(pi) = -{diA^) + -L{diA,{p,~Aj) + {p,^A,)d,Aj), (45) 
at at 2m 

whence 

(Pi h 

^j^^i^i) = i^^) + -^ijkiidjS - Aj)Hk) + —EijkidjHk). (46) 

In the case where p{t, x) is sufficiently localized compared to the characteristic scale of variations 
of the electromagnetic fields the angle brackets can be carried through the electromagnetic fields 
to obtain a closed system of evolutionary equations on the average position. They are the Newton 
equations with the stochastic correction. 

Now we return to the interpretation of the gauge fields A^. Recall that the expectation value 
of noise in the Langevin equation ( |3^ is called a systematic drift. In our case the systematic drift 
is equal to 

Mt,x):=diS{t,x)-Ai{t,x), (47) 



where we set ?ti = 1. Therefore to satisfy the second equation in (42) we have to take 



Ao-dtS = ^{f + hd,f). (48) 



Such an interpretation for the Langevin equation with a non-conservative force was proposed in [B4|. 
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The system of equations (^), (|4^ ) with respect to A^{t,x) and S{t,x) obviously admits a solution 
for any systematic drift fi{t,x). The fields A^{t,x) and S{t,x) are not uniquely defined by these 
equations and the arbitrariness in their definition is equivalent to the arbitrariness of a gauge. 
In particular, in the "unitary" gauge S{t, x) = the gauge fields Ai are the components of the 
systematic drift with an opposite sign. 



The Newton equations ( ^6] ) for the average position of the particle in the representation (47), 
become 



;^M = {/•). ;^W = {0. + /'8,)/'> + f 



M = {/•). i:^{':')^{(B, + f'dj)r) + -{Af). (49) 



For example, if f^{t,x) is the velocity field of an incompressible viscous fluid with the specific 



pressure p and kinematic viscosity i/ (see, e.g., |3^) then the second equation in (49) is replaced 
by 

^{xi) = -{d,p) + {v + h/2){Afi), (50) 

i.e., the acceleration of mean position of the Brownian particle is the same as for the particle which 
is not influenced by stochastic forces but entrained by a fluid with a higher viscosity. 

To gain a better physical insight into the stochastically deformed model of a nonrelativistic 
particle we construct the functional integral representation (^) of the transition probability. The 
Weyl-symbol of the Hamiltonian with the momenta pj + djS arising in formula (^) is 

Hit, X, ip) = -^ [V + 2ip\diS - Ai) + [diS - Aif] + A^. (51) 

Substituting this expression into (^) and integrating over momenta we arrive at 

Te(t,t') 

{t'-dr 
-i / dT [^±^ + {Ai - diS)x' - [A'' + drS)\ >, (52) 

where the functions Afj_{t, x) and S{t, x) obey the quantum Hamilton- Jacobi equation (|4^ ) and are 
taken at the point {t,x) = {t,x{t)). Since the expression in the exponent is the classical action 
modulo boundary terms the main contribution to the transition probability is made by the paths 
approximating a classical trajectory that is in agreement with our general considerations. In the 
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representation p^), (|^) the transition probability ( |5^ ) reduces to the well-known result 



t'~dT > 

V. (53) 



re{t,t') 

t'~dT 



So we have obtained the transition probability by the use of the Weyl-synibol of the evolution 
operator. The transition probability in terms of the pf^-symbol (31) of the evolution operator is 



constructed along the same lines and looks like (p^) without the stochastic correction to the action 
and the midpoint prescription. 

In conclusion of this section we briefly discuss how different generalizations of the Fokker-Planck 



equation of the form (^2D can be constructed in the developed framework. 



2. Fokker-Planck equation with nontrivial diffusion tensor 

At first we consider stochastic deformation of the model (35) on a curved space background 



with the inverse metric g'^K From the general considerations we know that the inverse metric is 
proportional to the covariance (the second cumulant) of noise in the Langevin equation (|32|), while 
the gauge fields Ai are related to expectation values of the noise. Besides, to provide a convergence 
of the integral (^) over momenta we have to require that the eigenvalues of the matrix g^^ are 
nonnegative. 

In constructing the Hamiltonian by its symbol we follow the simplest prescription [^] based on 
the use of the exponential map from the tangent bundle to the configuration space generated by 



the Levi-Civita connection (for more sophisticated methods see, e.g., [37, ^]). In other words we 
recognize the momenta pi as the derivatives at the origin of the Riemann normal coordinates and 
the vectors of the linear space V as the scalar functions on the manifold. Neglecting for a while 
the gauge fields A^ we have for the Hamiltonian (p^) in the normal coordinates 

^0 = ^ {PiPj 9'^ + '^Pid'^Pj + g'^PiPj ) • (54) 

Making use of relations between the derivatives of the metric taken at the origin of the normal 
coordinates and the Riemannian tensor [B^, 4C] we arrive at 



Ho = — (v'--], (55) 

2m V 12 ' ^ ^ 



We use the sign convention for the curvature as in M, M. 
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where R is the scalar curvature and V,- is the Levi-Civita connection. Notice that the covariant 



pg-ordering leads to the coefficient 1/3 at the curvature as in [40|, the covariant gp-ordering gives 
rise to the vanishing coefficient (minimal coupling) and 

corresponds to (|55| ) with the coefficient 1/6 at the curvature (conformal coupling) |^l|. All of these 
prescriptions result effectively in a variation of the coefficient at the scalar curvature. 

Thus the Hamiltonian for a particle interacting with the gauge fields A^ on a curved space 
background looks like 



H 



{hVi + Aif^—- 
2m 24771 



+ ^0. 



(56) 



Substituting the phase definition (^) into the Schrodinger equations following from the action (41) 



we obtain |17| 



g 



"'dt{g'"p) 



-V 



Zm m 



a.5 - ^„ + <*^ = -«, ln»V^ - Av.(a.S - A) + 1^. 



(57) 



where p and S are assumed to be the scalar functions and g := det gij. Since the Fokker-Planck 
equation is independent of the scalar curvature term the averages of observables depending on x 
only and their evolution do not depend on the concrete covariant ordering scheme. 

Rewriting equations ( |57[ ) in terms of the density g^'^p one can see that the inverse metric is 
proportional to the diffusion matrix. Besides from the Heisenberg equations we find a systematic 
drift 



mj^{x') = {g'HdjS - A,) + ^g-'/^d,{g^/'g'^)). 



(58) 



Despite the second term looks noncovariant it behaves under general coordinate transformations 
like the ffi'st term under the assumption that p tends to zero at spatial infinity. The transition 
probability @) takes the form (cf. ^ H) 

re{t,t') Te[t,t') 



Git',x';t,x)=[ n dMr) U ^, 



t'-dr 



exp 



h 



dT 



PiX H 

m 



-^ + ^p, U\d,S - A,) + \g-'l^d,{g'l\^^) 



(59) 



It is covariant under general coordinate transformations. To prove this by making a change of vari- 
ables one should take care about the so-called extraterms (see, e.g., [@^]) which cancel derivatives 
of the Jacobian matrices resulting from the noncovariant expression in the exponent. 
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3. Klein-Kramers equation 

Now we turn to another generalization of the Fokker-Planck equation (p4 ) , namely, to the Klein- 
Kramers equation. This equation arises in studying systems of second order stochastic differential 
equations. Therefore, firstly, introducing additional variables we reduce such a system to the 



system of the first order equations of the form (32) and then apply the developed formalism. 



Let us consider a system of stochastic equations 

x' = v\ v' = f{t,x,v)+u\ (60) 

where i^* is a Gaussian white noise. In that case the points of the configuration space are the pairs 
{x,v). We denote by ip,ir) canonically conjugate variables to {x,v) in the phase space 

{x\pj} = 6i, {v\7r,} = 6i. (61) 

From the representation (^Tj) it is not difficult to see that this stochastic system is described by 
the Hamiltonian 

H=^^^^^ + v% + A^ (62) 



where Af^{t,x,v) are the gauge fields. The Hamiltonian ( |62D is merely a general expression" at 
most quadratic in momenta ttj and linear in momenta pi. By introducing the phase S{t,x,v) the 



Schrodinger equations corresponding to the Hamiltonian (32) can be cast into the form 

h 



dtp = -d] 



^d^p+{drS-A,)p 



dlivip), 

(63) 



dtS -A^+ ^ - ^ '^ + v'd^S = --dlidfS - A,). 

The first equation is the Klein-Kramers equation, while the second equation is an analog of the 
Hamilton- Jacobi equation or, from an utilitarian viewpoint, is the definition of Aq. These equations 
are Lagrangian with the Hamiltonian action of the form (|4l]). Besides they are invariant with 
respect to the gauge transformations 

S{t, X, v) -^ S{t, X, v) + ^(t, X, v), Ai{t, X, v) -^ Ai{t, x, v) + d^^{t, x, v), 

Ao{t,x,v) -^ Ao{t,x,v)+i{t,x,v) +v^dfi{t,x,v). (64) 



A specific form of the coefficients at pi is fixed by the classical equations of motion for x^{t) which should coincide 
with @. 
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The covariant derivatives respecting these gauge transformations read as fohows 



'Po =Po + v'pi - Ao, 



Vi = TTi - Ai. 



(65) 



As far as an interpretation of the gauge fields Ai is concerned imposing the "unitary" gauge S = 
we see that these fields are equal to the components of the systematic force /j with an opposite 



sign. A path-integral representation of the transition probability (31) reads as 

d'^p(r)(i'*7r(r) 

rG(t,t') -rg[t,t') 



Git',x';t,x)= f n {dUr)dMr)) U ^^^ 

rg(t,t') rgft.t') ^ ' 



t'-dr 



exp 



h 



dT 



IT 



Pix' + TTii,' + i(-Y + iT^'{diS- Ai) + ipiv' 



■ (66) 



For the Hamiltonians at most quadratic in variables, i.e., for the linear systematic force /* ]44], 
an evolution of the stochastic system can be easily found from the Heisenberg equations. Here we 
consider the simplest case 



r 



-jv 



(67) 



just to illustrate how the formalism works. For the force ( |67| ) we can choose 



,t;2 n 



S = -7y , Ai = ^ Ao = j'-- -7d, (68) 

whence in the Heisenberg representation it follows from the definition of the phase S that 



{ip\7r,{0) = -7{^\v,iO), {ip\pi{0) = 0. (69) 

A general solution of the Heisenberg equations looks like 

p{t) = p, x(t)=x + 7~ (j5t — vr) — 7~ (-0 — 7" p) sinh7t + 7" 7rcosh7i, 

7r(t) = 71" cosh 7t + (7-0 — 7~ p) sinh7t, v{t) = 7" p + (£> — 7" j5) cosh7t + 7" -n-sinh7f, 

(70) 

where all the operators at the RHS of these equations are taken at i = 0. Making use of equa- 
tions (68), ( |70| ) and the commutation relations one can obtain an evolution of the average of any 
observable. For instance, 



9 , , , mh , f m , 9, mh 
^^^^^ - d+i—{v^) 



-{v\t)) 
2^ ^ " 47 



d\e-^^K 



(71) 



2' ' 47 

Then assuming that 7 > and the equipartition law is fulfilled we determine the deformation 
parameter in equations ( |63|) 

2-ikT 



h- 



m 



(72) 
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It is interesting to note that the trace of the equipartition law in the Fokker-Planck equation (42) 
with the Hamiltonian (^) is the relation 

hm n^dt} = ^^[x. [.%H]] = ^d, (73) 

where x{t) are the position operators in the Heisenberg representation, T means the chronological 
ordering and hpp is the deformation parameter in the Fokker-Planck equation. To reproduce the 
equipartition law one should formally put dt = 2^~^. 

B. Relativistic particle 

In this subsection we consider three relativistic models which under stochastic deformation give 
rise to relativistic generalizations of the Fokker-Planck and Klein-Kramers equations, and to the 
equation describing a massive charged particle influenced by external systematic and stochastic 
forces with the radiation reaction taking into account. In other words in the last case a dissipation 
is described by the Lorentz-Dirac force |45, E^]. 



1. Relativistic Fokker-Planck equation 

In the previous subsection we saw that the stochastically deformed model of a nonrelativistic 
particle results in the Fokker-Planck equation. Therefore it is reasonable to expect that stochastic 
deformation of the model of a relativistic particle gives some relativistic generalization of the 
Fokker-Planck equation. 

The Hamiltonian action of an interacting relativistic particle has the form 

Sh[x,p,X] = Jdrlp^x^^ + ^(P2 _ ^2)]^ (74) 

where 7^^ '■= Pfj. — A^ and A^ are the gauge fields. The dynamics of the model (^) are governed 
by one first class constraint. In the proper time parameterization A = m~^ the evolution of the 
position of the particle obeys the equation 

mi;^ = ^^ - p^. (75) 

Thus by analogy with a nonrelativistic model one should expect that the corresponding stochastic 
equations look in the unitary gauge like 

mx^" = A^' + ..., (76) 
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where dots denote terms vanishing at h = 0, i.e., in the classical limit A^ plays the role of the 
4-momentum of the particle. The problem that one encounters by naive introduction the noise in 



the RHS of Eq. (76) is to preserve the reparameterization invariance. In the proper time gauge, 
for example, we must guarantee the fulfillment of i;^ = 1. The stochastic deformation procedure 
allows us to obtain the Fokker-Planck type equation associated with ( [76| ) respecting the gauge 
invariance. 

We shall deform the model ( [7^ ) in the gauge of a laboratory time 

r = x°. (77) 



A thorough description of quantization of a relativistic particle in the gauge (77) is presented in 
ipTII and we just trace some basic steps of this procedure which are necessary for us. Firstly, we 
solve the mass-shell constraint with respect to the energy^ 



Vo - ^m2 + (Pi)' = 0. 



(78) 



Then we naturally realize the Heisenberg-Weyl algebra of operators in the linear space V of two- 
component real vectors 



^ 



tp{x) 



GK 



rt^ .- 



P^.■■-- 



-hd„ 











-hd,. 



(79) 



xf" 
x^' 

Define the dual linear space V* as the space of real linear functionals acting on elements of V by 
the rule 

'"o 1 

1 



($1*) := 



dx$^ 



* 



dx[c^^' -I- ip'i/j], 



(80) 



where 



$ 



(fix) 



(81) 



3D. 



The pure states are as usual the projectors of the form^*^ 

The above construction allows us to realize the constraint ( |78| ) without going to pseudo- 
differential operators: 



T:=Vo 



m2 + {Ai - hdif 



(82) 



Here we work in the particle sector only. Antiparticle's states in the sense of [M7| are not included. 

As in Section O one can construct elements of the linear spaces V and V* formally acting by the operators p^ on 

the doubly degenerate vacuum of the operators x'^ . 
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The dynamics of a pure state of the stochastic system are governed in the gauge (|77|) by the two 
Schrodinger equations 

f|^) = ($|f = 0. (83) 

In quantum mechanics they are simply the Klein-Gordon equation in the first order formalism |48]. 
In our case these equations are equivalent to 

{V^ - m'^)'iP = 0, ilj' = Voip, {V+^ - m^)ip = 0, ^' = V^^, (84) 

where the cross denotes a formal conjugation. When the gauge fields vanish these equations are 
the Klein-Gordon equations for tachyons^^ . On the solutions of the Schrodinger equations ( p3|) the 
inner product (pO) reduces to 

($1$) = - j dx[<fVoi^ + i^V+^]. (85) 

The probability density function corresponding to the state |$)(^| is proportional to the integrand 

The Heisenberg equations associated with the Schrodinger equations (p^) allow to obtain the 
evolution of an average of any physical observable and, in particular, to establish a stochastic 
interpretation of the gauge fields A^. However we shall act in a different way to derive a direct 
stochastic interpretation of Eqs. (|83) in terms of the Langevin equation of the form (|3^). To this 
aim we, firstly, represent the fields (p{x) and 'ip{x) as 

ip{x) = e^^(^), ^(x) = p(x)e-^^(^), (86) 



substitute them into Eqs. (pi) and fix the unitary gauge 5 = 0. Then we have 
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Qf^ 






0, A^ -m^ = hd^Af,. (87) 



The first equation is the conservation law of the probability density current 

JM = ^(^M + f5M)p> (88) 

while the second equation is the deformed mass-shell condition or the mere definition of A^{x). For 
given Ai(x) we take the solution of the second equation in ( [87| ) which is regular in h and possesses 
the following classical limit 

^oU=o = Vm^ + (A?- (89) 



^^ For the interrelation between relativistic random walk models and relativistic wave equations see, e.g., |[43, pU J . 
^^ Notice that for the constant fields A^ Eq. mm reduces to the well-known relativistic diffusion equation J18[ . 
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The solution defined in this way is unique. Consequently, if one has jfj,{x) at the initial moment 
then one can find from Eqs. (|8q ) the function p{x) and its first derivative with respect to time at 
this initial moment and solve the Cauchy problem. 

For the potential and stationary relativistic flow with a constant specific relativistic enthalpy 



|35| entraining the Brownian particle we can choose Ai = diU{x), where [/ is a Lorentz scalar. 



Then the stationary probability distribution function following from ( p^ ) looks like 

f = Vl + (Ai)2/m2 exp{-2U/h)/Z = jexp{-2U/h)/Z, Z = const, (90) 

where 7 is the Lorentz factor. The power of the 7-factor is uniquely specified by the requirement 
that j^ is the zeroth component of the 4-vector. 

It is easy to see from Eq. (p8|) that the derived equations (84) lead to negative probabilities 



when the probability density function changes rapidly on the spatial (time) scales less than or 
comparable with the "Compton wave length" hm~^. This is in a perfect analogy with the well- 
known property of the Klein-Gordon equation in relativistic quantum mechanics [51|. Therefore 



Eqs. (|84| ) have a correct stochastic interpretation only for the fields with a characteristic scale of 
variations much larger than hm^^. 

Despite the above remark we consider an evolution of the probability density function localized 
in one point x* = y' at the initial moment x^ = y^: 

f{y',x') = 6''-\x'-y'). (91) 

The solution of the Cauchy problem to the first equation in (^) is obtained via a convolution of 
the retarded Green function with the initial distribution 

l + ^apo]6''{x-y), (92) 

where the function a(x) is the regular in h solution of the equation 

^0 + ^do) a(x) = 1. (93) 

Thus, formally, the transition probability is 

G(x, y) = (ao - ipo) T^^^T^ (1 + l-Po) 5% - y)- (94) 



The transition probability (94) does not possess the defining property of the Markov process 
and it can not be straightforwardly represented by a path-integral. To overcome this difficulty we 
represent the transition probability (|9^) as the sum of the transition probabilities possessing the 
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Markov property with respect to some new variable (fictitious time). The stated procedure is an 
analog of what is called in relativistic quantum mechanics the Fock's fifth parameter (proper time) 



formalism [E^, 53] 



Under some technical assumptions on the transition probability (^) the following equality for 
the operator G with the kernel (p3) holds 



1/1 \ / 1 \ -1 



G = h-^ j dre-i^, H:=(l + ^apo) Qp^ - p'^A^ (aq - ^po) • (95) 



This formula is understood in the distributional sense. It is derived by an application of the descent 
method (see, e.g., [|^, |2^) to the equation 

- hdripir, x) = Hi;{T, x), ^'(O, x) = 6'^{x - y), (96) 



and is valid when the integral (|95D converges. It is not difficult to prove the convergence of this 
integral in the case of the constant fields Ai(x). In a general case the integral ( [95| ) can be defined 
at least perturbatively. 

The evolution operator exp [—h^^rH] respects the Markov property and conserves the proba- 
bility density function normalization. Its kernel in the x-representation can be interpreted as the 
probability of particle's arrival to the point x of the space-time from the point y with the value of 
the fictitious time r. The sum over all such particles with different fictitious times gives the total 
observable transition probability (p3). 

Now we represent the kernel of the evolution operator exp [—h^^rH] by a path-integral along 
lines expounded in Section || and introduce the corresponding Langevin equation (|3^). It seems 
impossible to find an explicit simple expression for the probability density function of noise for 
arbitrary fields Ai{x) but some properties of this distribution are easily derived. 

Firstly, in the classical limit /i ^ the Hamiltonian (^) reduces to 

H = -p^A^/Ao, (97) 

where ^o is specified by Eq. (p^) . Consequently, the probability density function of noise reads 



F{iy, r, x) -^ 5(i/° - l)^''-^ u' - a' /^/m^ + {Ai)^) . (98) 

We see that in the classical limit the fictitious time turns into the physical time x^ of the lab 
frame. Besides the above formula gives the interpretation to the fields Aj(x), which confirms our 
expectations. 
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Secondly, in the nonrelativistic limit Aq ^ m and m ^\\pq\\ <^ 1, where po is assumed to act on 
the probability density function, the Hamiltonian (P5|) is rewritten as 



H = Jlil.i!l^, (99) 

2m m 

whence the standard Gaussian distribution for the noise follows. As in the classical limit the 
fictitious time r is the physical time x^. 



Thus Eqs. ( |84D generalize the Fokker-Planck equation (42) in the same sense as the Klein- 
Gordon equation generalizes the Schrodinger equation. In the classical limit these equations de- 
scribe the same physical system and in the nonrelativistic limit Eqs. (p4|) reduce to the Fokker- 



Planck equation (42). However it is worthy to note that the above interpretation in terms of the 
fictitious time has some weaknesses. Namely, since we start from the probability density function 
(|9l| ) localized in one point the obtained equations are ill-defined. This inevitably leads to negative 
probabilities although they disappear in the classical and nonrelativistic limits, when hm~^ tends 
to zero. Moreover the negative probabilities appear to result in the probability density function of 
the noise acting on the particle arrived to the point x with the fictitious time r being not equal 
to zero in the superluminal region. But, of course, the sum over all the particles with different 
T gives rise to the transition probability with the support bounded by the light cone. A detailed 
investigation of these peculiarities will be given elsewhere. 

2. Relativistic Klein-Kramers equation 

To formulate a relativistic generalization of the Klein-Kramers equation we start from a classical 
relativistic system with the Hamiltonian action of the form (cf. (p3)) 



Sh[x,p, y, vr. A, /i] = / dT{p^x^' + vr^y^ - ATi - 1JT2), 



Ti:=^gt^-{Tr^ + A^)iTr,+A,)+y''p^ + (j), T2 := y^ - m^ (100) 

where dots denote the derivatives with respect to the parameter r, m is a mass of the particle, 
A^{x,y) and (f){x,y) are gauge fields, g^'^{x,y) is a symmetric Lorentz tensor orthogonal to y^. 
In other words the dynamics of the system are governed by two first class constraints. The first 
constraint generates reparameterizations of the world line of the particle, the second constraint 
being just the mass-shell condition. These two constraints are in the Abelian involution by virtue 
of the property of the tensor g^^'^ . 
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A part of the equations of motion in what we are interested in reads as follows 

i;'^ = A2/^ r = ^9'"'{A, + ^,), y^ = m\ (101) 

By analogy with a nonrelativistic case we expect that in the classical limit ^ ^ and in the unitary 
gauge 5 = the stochastically deformed system ( |100D will describe a relativistic particle with the 
equations of motion 

y>^ = m-^g^^'^ix, y)A,{x, y)V^, m^ = yf'. (102) 



Hence in the unitary gauge and in the classical limit m^^g^^'^ Ay is expected to be a systematic force 
acting on the particle. As in the previously considered model a naive introduction of the noise to 
the RHS of the first equation in ( |102D spoils the reparameterization invariance. To preserve it we 
have to demand from the noise correlators 

y^^{T){u^{T)v^,{n) . . . u^^Tn)) = 0, (103) 



where Vjj.{t) is the noise. The identities ( |103D can be viewed as the Ward-Takahashi identities for 
this model (see a detailed discussion in |^, ^). By the use of the stochastic deformation procedure 
we shall find below the general form of the Fokker-Planck equation associated with ( 102D respecting 



the reparameterization invariance in the case of the Gaussian white noise. A generalization of the 
obtained equation to the case of an arbitrary noise distribution law will be obvious. 



To deform the model ( 100 ) we naturally realize the Heisenberg-Weyl algebra in the linear space 



V of smooth functions ijj{x,y): 

P^ = -hd^^, r = y^ ^^. = -hdy- (104) 



Then the constraints ( |100| ) are the Weyl-ordered operators 

fi = hu,X9^'' + 2U^,g^''Uu + g^''U^U,) + y^'p^, + <P, f2 = y^-m\ (105) 

o 

where H^ = tt^ + A^. These constraints are in the Abelian involution as well^^. The physical pure 
states of the system are singled out by the conditions 

y^^fi\ilj) = {^\yQ^fi = 0, fblV) = {v\f2 = 0, (106) 



^^ There are simpler forms of the first constraint which preserve the classical algebra of constraints. For example, 
one can choose 



Nevertheless this arbitrariness is unimportant and affects the definition of (j) only (110) 
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where the inner product is defined in a standard way and we resolve the first constraint with respect 
to pq. Equations (|106|) are invariant under the gauge transformations 



Besides the transformation of the gauge fields A^ of the form 

Ati{x,y) -^ Af,{x,y) +y^C{x,y), 



y^'dU. 



(107) 



(108) 



leaves Eqs. ( |106D unchanged. Therefore the fields A^ have d — 1 independent components and in 
the massive case one can make them orthogonal to y^. Introducing the phase 



^(x,y) =ei^(^'^), ^(x,y) =yV^,y)e-i^(^'^), 



(109) 



we reduce the first pair of equations (the Schrodinger equations) in formula ( p. 061) to the system 

'h 



QV 



-g^'dy-p-g^^idis + A^Yp 



dl{y^~p)=^, 



\g^'{dlS + A,){dyS + A,) + y^d^S + 4> = -^dlig^^idlS + A,)] - %dl,g^^^ 



(110) 



The first equation is the conservation law of the probability density current whose time component 
is the probability density function. The second equation is the definition of (p. 
The second pair of equations in ( |10£| ) is taken into account by the equalities 

p = 6{y^-m^)p{x,y'), S = ln6{y^ - m^) + S{x,y'). (Ill) 

On substituting these equalities to the first equation in ( |110| ) we obtain a relativistic generalization 
of the Klein-Kramers equation |2l|, ||, |2|, ||, ||, |2|, ||]i^ 



Of 



^05|p_(5|5 + ^,)p 



^^ ( V 



0, 



(112) 



where the positive root of the mass-shell condition is chosen: y^ = ^Jvrfi -|- {yif' ■ The logarithm 
of the (5-function in (|111| ) does not actually contribute to Eqs. (|112| ) because of the orthogonality 
of the tensor g^^ to y^. Now it is not difficult to construct a path- integral representation of the 
transition probability and the corresponding Langevin equations (32). For brevity we give the 
Langevin equations only 



x^ 



(^^(r))=0, (^y\r)v\T')) = h^^b(T - t'\ 



(113) 



^"^ The position of the tensor <;'■' in Eq. (112) between two derivatives 9f is uniquely determined by the requirement 
that the evolution keeps the mass-shell condition. It is that order which results in the Jiittner stationary distribution 
derived from the relativistic Klein-Kramers equation in the papers vA. 
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where z^* is a Gaussian white noise and we assume that g^^{x,y) is positively definite. Recall that 
the stochastic equations (|113D are understood in the Ito sense. Thus we confirm our expectations 



that equation ( |112| ) or ( |113D describes a relativistic particle in the gauge of a laboratory time (77) 
under the influence of noise. 

As far as the tensor g^^'^ is concerned its physical meaning is comprehended from the Langevin 
equations ( |113| ). We remark only two possible choices 

5^ = ^-ry^^ 9!^'' = -{ri',-^W-^), n' = l. (114) 

The first tensor implies that the effect of the stochastic force on a relativistic particle is isotropic 
in a momentary comoving frame. The second tensor corresponds to the isotropic influence of the 
stochastic force in the frame distinguished by the vector n^. It is interesting to note that in the first 
case, which is minimal from the mathematical point of view since new objects are not introduced, 
the tensor g^^ can be interpreted as the inverse metric in the momentum space with the coordinates 
2/*. Keeping in mind that 

(det9i^y/2 = m-^y0, (115) 

we can rewrite Eq. ( |112| ) in a covariant way. The metric g]- is a metric of a constant positive 
curvature and, consequently, it is conformally flat. The second tensor g2 becomes Euclidean 
provided that 71^ = 5^. 

To conclude this example we briefly discuss a representation of the transition probability in 



terms of the proper time. The main observation allowing to formulate Eqs. (|113D in the proper 
time gauge is the following. The retarded Green function of the twisted forward Schrodinger 
equation (|106|), that is the transition probability, is the kernel of 



oo 



^i = -y^^^— = [ dT^ei^^^^^-\ (116) 

where (p := y^ ip and the last equality is proved by the descent method as it was done in consid- 
ering a relativistic generalization of the Fokker-Planck equation. The operator in the exponent is 
proportional to the operator acting on p in the first equation in formula ( [LlOj ). 
The kernel 

(x,y|exp -^(ffiif^^ \x',y') (117) 

is the probability density of particle's arrival to the point (x,y) from the point {x',y') with the 
proper time r measured in the momentary comoving frame. The factor y^m^^ in (|116| ) is caused 
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by the passage from the momentary comoving frame to the lab frame. The transition probability 
possesses the Markov property with respect to the proper time r. It conserves the normalization 
of the probability density function, i.e., the probability to find a particle with the proper time r in 
some point of the space-time with a certain momentum is equal to unity. Since the generator of 
evolution commutes with the operator T2 the transition probability ( |117D respects the mass-shell 
condition. 



The Langevin equations associated with (117) read as 



(118) 

whence we infer that r is indeed the proper time. As before the noise i'^ is a Gaussian white 
noise. The representation ( |116| ) of the transition probability is especially useful when a relativistic 
particle has a finite lifetime. In that case the integral (|116| ) over r is cut on the upper limit by the 
lifetime of the particle provided, of course, that the particle was created at the initial moment. 

3. Stochastic Lorentz-Dirac equation 

A relativistic generalization of the Klein-Kramers equation regarded in the preceding example 
does not apply to a charged particle with a radiation reaction taken into account. If we apply the 
above scheme to the so-called Landau-Lifshitz equation |^], which is obtained from the Lorentz- 
Dirac equation by the reduction of order procedure, then the vanishing external systematic force 
entails a zero dissipation force in the resulting stochastic equation. That is why we should stochas- 
tically deform the Lorentz-Dirac equation itself. The procedure of stochastic deformation is very 
similar to what we have considered in the previous example. Therefore we mention the crucial 
points only. 

We start from the Hamiltonian action of the form 

Sh[x, X, y, y, w, tv, A, /U, x] = / dT{x^x^' + y^y^" + w^w^" - XTi - fiT2 - xTs), 

Ti := ^5^"ir^VF. - loV^M + «^% + y% + </'' T2:=y^-1, T^.= yw, (119) 

where W^ := w^ -|- A^ is a covariant momentum, A^{x^y^w) and (j){x,y,w) are gauge fields and 
g^'^{x,y,'w) is a symmetric Lorentz tensor orthogonal to y^. For simplicity we assume that the 
particle has a unit mass. The constraints Ti, T2 and T^ are of the first class with the algebra 

{Ti,T2} = -2T3, {Ti,r3} = w^T2. (120) 
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A relevant part of the equations of motion resulting from the action ( [119[ ) reads as 

x'' = Xy^', y'' = Xw'', w'' = Xig'^'W^ - w^y''), y^ = l, yw = 0. (121) 

These equations suggest that in the classical limit and in the unitary gauge the stochastically 
deformed model ( |lig| ) describes a relativistic particle obeying the equations of motion 

x^' + x'^x'' = g^"'A^, x^ = y^, y^' = w'', (122) 

where the derivatives are taken with respect to the proper time. The LHS of the first equation is 
proportional to the Lorentz-Dirac force and we refer to this equation as the Lorentz-Dirac equation. 
The mass term is contained in the RHS of this equation. 

The Heisenberg-Weyl algebra is naturally realized in the linear space V of smooth functions 

x'' = x'', 3:^ = -hd;, r = y'', y,^ = -fidy, w^' = w'', i!^ = -hd;:. (123) 

The constraints turn into the appropriate Weyl-ordered operators. Their algebra coincides with 
the classical algebra of constraints ( |120| ). The physical states of the stochastic system are specified 
by 

y^'tl\^l;) = {^\y^'fi = 0, fslV^) = ((^ifs = 0, fslV) = ((^iTg = 0, (124) 



where the standard inner product is understood. The first pair of equations in formula ( |124| ) is 
just the forward and backward Schrodinger equations. 

Making a substitution of the form (109) into the first equation in (124) we arrive at^^ 






l9'"'d::p - {g'^'id^S + A,) - w^yf^) ~p 



dl{w^p)-dl{y^^~p) = Q, 



^5''"(5;r5 + A^){d'^S + A,) - w^y^'d'I^S + w'^dyS + y^'^^^S - yw + cj) = 

- ^9-b^'^(9-5 + A,)] - ^d;:,g^r (125) 

As usual the first equation is the conservation law of the probability density current. The prob- 
ability density function is the time component of this current. The rest of equations defining the 
physical state results in 

p = 5{y^ - l)5{yw)p{x, y\w'), 5 = In [6iy^ - l)5{yw)] + S{x, y\ w'). (126) 



The proper time Liouville equation for the Lorentz-Dirac equation can be found, for example, in |58 
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Under the above conditions the first equation in (|125|) can be reduced after a httle algebra to 



df 



g'^df -{g'^{dfS + Aj 



w'y^] 



vl 



df I w 



^M ( y'i 



Vo 



0, 



(127) 



where we have assigned y'^ = ^1 + (yj)^ and w'^ = y^ yiWi. Recall that the probability density 
function is y^ p. Now a path-integral representation of the transition probability associated with 
Eq. (127) is straightforward. The Langevin equations (32) with the Ito prescription become 



w 



x'^ 



w 



9 



2y 



h. 



y^{dJS + A,)-w^^^+-dJ 



1^1 



+ v' 



-9 



•*j 



(128) 






where ly^ is a Gaussian white noise and g^-'{x, y, w) is positively definite. As regards explicit forms 
of the tensor g"^^ see the discussion in the previous example^^. The obtained stochastic equations 
are equivalent to the Lorentz-Dirac equation ( |122| ) with a random force which is rewritten in the 
gauge of a laboratory time, that is equation ( |128| ) is the stochastic Lorentz-Dirac equation. 

The proper time representation of the transition probability from the point {x',y',w') to the 
point (x, y, w) takes the form 



{x,y,w\ I drh y exp —(pTi(p 



\x ,y ,w} 



(129) 



where (p := y^ tp and (pTi(p~^ is the operator acting on p in the first equation of ( [1251) . Besides, 
the initial state should satisfy the mass-shell condition and its differential consequence. Therefore 
the proper time stochastic Lorentz-Dirac equation is the following 



x^ = y^ 



ij^ =w 



^ w^^ = g^''{d:S + A,) - w^y^ + ^d^ g^^^ + z.^ 



{u^{t)) = 0, {v^^{t)v-{t')) = hg^-5{T - r'), 

where v^^ is a Gaussian white noise. 



(130) 



C. Free relativistic field models 

In this subsection we consider the stochastic deformation of two free models: a relativistic real 
scalar field and an electromagnetic field. We regard the scalar field model as a prototype of any 



In the series of papers |M, pA the stochastic Lorentz-Dirac equation was derived from quantum electrodynamics 
in the paradigm of the Feynman- Vernon influence functional. There was obtained an exp licit expre ssion for the 
tensor 3^" that is a nonlocal in x^ operator orthogonal to j/p. A generalization of Eqs. (125) and (127) to this 
case is straightforward. 
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field model. Besides, it is the simplest model of a one-dimensional crystal in the continuum limit, 
where the scalar field describes the displacements of points of a medium from their equilibrium 
positions. The stochastically deformed model of Maxwell's fields allows us to touch the problem 
of electromagnetic fluctuations [^] bringing it in accordance with the general scheme advocated 
in this paper. 



1. Scalar field 

Consider a real scalar field 0(x) on the Minkowski space M^''^^-^ with the dynamics governed by 
the action 



5M = iy"d'^x(a^</.9^0-mV), 
where m is a mass. The corresponding Hamiltonian action has the form 



(131) 



d'^x 



where the dot denotes the derivative with respect to time. The Hamilton equations are 

(j) = IT, TT = A(p — m (j). 

The Noether theorem gives the conserved 4-momentum 

Vo = H = - / dx(7r2 + d^(j)di(p + 771^02), ^i = / d:!C7rdi(/), 



(132) 



(133) 



(134) 



which is the generator of translations. 

Now we are going to define the stochastic deformation. The canonical coordinates (j){x.) and 
canonical momenta vr(x) turn into the generators of the Heisenberg-Weyl algebra 



(x),7r(y)]=M(x-y) 



(135) 



Since we deform the linear model its Heisenberg equations coincide with the Hamilton equations. 
The Heisenberg equations can be formally solved by means of the expansion in terms of the complete 
set of solutions of the classical equations of motion in the form of plane waves. The field operator 
and its canonical conjugate read as 

dp 1 



n{x) 



(27r)^-i ^ 

dp 
(2^0^^ 



a(p) cos(p^x^) + 6(p) sin(p^x^) 
b{p) cos{p^x^) - a(p) sin(p^x^) 



(136) 
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where po := \/vn^ + p^. As it follows from ( |135| ) the operators a(p) and 6(p) obey the commutation 
relations 



[a(p),6(k)] = ;i(2^r-^5(p-k). (137) 



The Weyl-ordered operators of the generators of translations ( |134| ) look like 



^^=/ (2#ryr(p) + ^'(p)]- (138) 



The solution ( |136| ) also implies that 

\W)Mv)\ = n I ^P^I^M^-^)'] = -2/.sign(x0 - y^)G{x - y), (139) 

where G{x — y) is the symmetric Green function (see, e.g., ]59|]) 

{U + m'^)G{x) = 5^{x), G{x-y) = G{y-x). (140) 

As in the case of quantum field theory two field operators separated by a space-like interval com- 
mute. 

The key ingredient of stochastic field theory is the propagator, which we define in a standard 
way 

G{x,y):={T{4>{x)4>{y)}), (141) 

where T is a chronological ordering. Any other correlator can be expressed in terms of the propa- 
gator by means of the Wick theorem and the commutation relation (139). 



Certainly, an explicit form of the propagator depends on the state entering in its definition. 
Nevertheless the propagator possesses the properties following from the rule of differentiation of 
the chronological ordering and the Heisenberg equations 

(n^ + m2)G(x,y) = -M^(x-y), G{x,y) = G{y,x), (142) 

regardless of the state which does not depend on time explicitly. Therefore the propagator can be 
represented as 

G{x,y) = -hG{x-y) + ..., (143) 

where dots denote some solution of the Klein-Gordon equation. 
Moreover if the state p is translation invariant 

[/5,^/.]=0, (144) 
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then the obvious identity 

Sp([P^,pT{0(x)<^(y)}]) =0 (145) 

entails a translational invariance of the propagator, viz. 

G{x,y) = G{x-y). (146) 

Analogously it can be proved that an invariance of the state with respect to the Lorentz transfor- 
mations results in a dependence of the propagator on the space-time interval only. 

Unlike quantum field theory there is no distinguished ground state in our case since the Hamil- 
tonian ( |138[ ) has an unbounded spectrum. That is why we define the state p in such a way that the 
propagator (|141D becomes minimal in view of the above considerations, that is the dotted terms 
in formula ( |143| ) vanish. 

Let us achieve this goal gradually. Firstly, suppose that p commutes with the Hamiltonian, i.e., 
the state is invariant with respect to translations in time. Then from the equality 

([a(p)6(k), i^]) = pMp)m) - A:o(a(p)a(k)) = 0, (147) 

and the same equality with p and k interchanged we find 

(a(p)a(k)) = (6(p)6(k)) = /(p, k)5(po - fco), /(p, k) = /(k, p). (148) 

Furthermore, from the equalities 



([a(p)a(k),i7]) = l(fco[a(p)6(k) +6(k)a(p)] +po[a(k)&(p) + 6(p)a(k)]) = 0, 
([6(p)6(k),iJ]) = -i(A:o[Kp)a(k) + a(k)6(p)] +po[Kk)a(p) +a(p)6(k)]) = 0, 



(149) 



we infer that 

^(a(p)6(k) + 6(k)a(p)) = ff(p, k)5(po - h), ff(p, k) = -^(k, p). (150) 

By the same way if we assume that the state p is also invariant with respect to the spatial 
translations then 

(a(p)a(k)) = (6(p)5(k)) = f{p)5{v - k), ^(a(p)6(k) + 5(k)a(p)) = 0. (151) 

Thus we have for the dotted terms in formula (|T4 



/(2^^^°^f^'^^"-^)']' ^'^'^ 
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where one should bear in mind that the propagator expressed in terms of the operators a and b 
looks like 



G{x,y) 



dpdk 



{ a(p) cosipf.xf') + b{p) sin(p^x^') o(k) cosik^yf") + 6(k) sin(A;^y'') ) , (153) 



(27r)2rf-2^^^H 

provided x^ > y^. The additional requirement of an invariance of the state with respect to the 
Lorentz group yields that the function f{p) contributes to the integral ([L52[) as a constant. Con- 



sequently, this integral is proportional to the so-called Hadamard function |59| 



G«(x-y):= 



dpcos[pf,{x-yy 



(154) 



{2ttY-^Po 

The Hadamard function does not vanish for the points x and y separated by a spatial interval. 
Hence if we want to obtain a causal theory we have to set the proportionality constant to zero that 
leaves us with the propagator proportional to the symmetric Green function. 

The above considerations prove that the propagator proportional to the symmetric Green func- 



tion is unique for the stochastically deformed model (131) under the causality condition and the 
requirement of the Poincare invariance of the state. Explicitly such a pure state p = \^){^\ can 
be specified as 



i^)('&i=ni^p)^^pi' 



(155) 



where 



[a\p) + 62(p)]|*p) = {%\[a\p) + b\p)] = 0. 



(156) 



In particular this implies 



P„|^) = ($|P„ = 0. 



(157) 



The propagator (141) naturally appears when one computes, for example, the transition prob- 
ability from the state |^i)(<I>i| to the state |*2)(*^2| 



($2|?7|*l)($l|?7-l|^2), 



(158) 



where U is the evolution operator over an infinite time and the states |^i)($i| and |^2)(^2| are 
obtained from the ground state ( |155D by the action of the operators a and b. The perturbation 
techniques is in a perfect analogy with such a procedure in quantum field theory. In our case 
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the so-called closed-time-path formalism (see, e.g., [^, 30 1 and references therein) seems to be 



especially useful for calculations. If the classical current is introduced into the model (131) then 
the solution ( |136 ) of the Heisenberg equations is the general solution of the homogeneous equation 
only. A particular solution of the inhomogeneous equation with current is given by a convolution 
of the retarded Green function with this classical current. In that case it is the convolution of the 
retarded Green function with current that is the average field over the ground state (ff55|). 



2. Electromagnetic field 

As usual we start from the classical Maxwell action-*^^ 

i2 na A A I ri^i A \2 



S[A] = -'^ f d^xF^.F^^'^ = ^ I d^x 



At - 2diAoA, + {d^AoY - {diAjY + diA,d,Ai 



(159) 



where F = dA is the strength tensor of the electromagnetic field. Making the Legendre transfor- 
mations on Ai only and introducing the canonical momenta 

dL 



TT,; 



dA^ 



diAo - Ai 



(160) 



we arrive at the Hamiltonian action 



Sh= f d^x JttM* - \ K' + {diA.f - d.AjdjA,] + A^diA . (161) 

As we see the time component of the gauge fields A^ is merely the Lagrange multiplier to the 
Gauss law constraint. The expression in square brackets with 1/2 is the density of the physical 
Hamilton function Hq. Thus we have a model with first class constraints. 

According to the standard BFV-quantization scheme for the first class constrained models we 
introduce the canonically conjugate ghost pairs (c, P) and (c, P) 



[c(x),P(y)] = [g(x),P(y)] = M(x - y), ghc = ghc = -ghP = -ghP = 1, 



(162) 



where square brackets denote graded commutators. Besides we introduce the canonical momentum 
ttq, ghvTo = 0, to the Lagrange multiplier ^o- Then we construct the BRST-charge 

iTTi + CTo), [n,n] = 2n'^ = Q, ghf^ = i. 







/ dx(c9j7 



(163) 



Physical observables and the Hamiltonian should commute with the BRST-charge. In our case the 
Hamiltonian Hq obviously satisfies this requirement. The gauge fixing fermion is taken in the form 

a , 



V' 



(ix 



PAq + P[ d,A, - -TTo 



ihtp 



-1, 



(164) 



We follow |l4] in the deformation procedure for the model of free electromagnetic field. 
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where a is an arbitrary constant. 

Therefore the gauge fixed Hamiltonian looks hke 

[d^l\ [^f + (9^4•)' - d^A,djA,j - Aodiir' - ttq (s^i^ + ^ttq^ - tP + d^cdiP\ , (165) 
and commutes with the BRST-charge as well. The corresponding Heisenberg equations read as 

Ai = diAo-TTi, Aq = diAi - ano, c = c, P = -AP, 

. (166) 

TTi = didjAj - AAi - diTTQ, ttq = -diiTi, c = Ac, P = -P. 

As expected the ghosts dynamics are split of the fields dynamics and we do not enlarge on them 
henceforth. Combining the Heisenberg equations it is not difficult to obtain 

aAo = ia-l)dm, aA, = {1 - a)di7To, Dvr, = 0. (167) 

In the Feynman gauge q = 1 the Heisenberg equations on the fields reduce to the wave equations 
and we can use all the results of the preceding example regarding the model of a scalar field. The 
physical states are those containing solely transverse photons. 

Hence in the Feynman gauge the Poincare-invariant causal propagator of the fields A^ becomes 

Gix-y) = {T{A^ix)A,iy)}) = hv^,Gix-y), nG{x)=6^ix), (168) 

where G{x) is the symmetric Green function, what coincides with the well-known result of [^] in 
the case of a transparent medium e = fi = 1. The correlators of the electromagnetic fields F^i, are 
straightforwardly obtained from (168). 



IV. CONCLUDING REMARKS 

There are, of course, many important open problems regarding the procedure advocated here. 
We mention a few of them only. Apart from the study both analytical and numerical of solutions 
of the derived relativistic equations we distinguish a more detailed investigation of the model with 
nonlinear phase space. In spite of the fact that we considered several models of this kind their a 
priori interpretation remains unclear. For instance, in studying the stochastically deformed model 
of a relativistic particle with Hamilton function quadratic in momenta and with a nonlinear phase 
space, we found after constructing a path-integral representation of the transition probability that 
the probability distribution function of the random force is non-Gaussian. It is desirable to obtain 
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a simple method to discover the noise distribution from the initial classical model. Furthermore we 
only touched the problem of stochastic deformation of field theories. The next step is to consider 
interacting models and their renormalization properties. It is also interesting to investigate the 
structure of the gauge transformations and corresponding functional gauge fields associated with 
the changing of the phase S in the stochastic field theory framework. Notice that the gauged 
formulation of stochastic mechanics developed in this paper allows us to represent it in terms of 
sections of vector bundles analogously to quantum mechanics. For the Schrodinger equations of 
the form (^) the pair of functions ip{x) and tp^x) describing a pure state of a stochastic system can 
be viewed as section (written in the light cone coordinates) of a vector bundle over the space-time 
with the typical fiber M^ and the structure group SO{l, 1). Then the states should be described 
by time- like or isotropic sections to guarantee the non- negativity of probability density functions. 
The gauge fields are connections on this vector bundle. A similar construction on the Whitney sum 
of tangent and cotangent bundles TM © T*M is called an almost generalized product structure 



(see, for introduction, ||6l[| ) that is a real analog of an almost generalized complex structure |62| 
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